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Dear  Tom: 


1  am  sending  you  an  annual  report  mostly  in  pictures.  I  chose  the  topics 
which  could  be  more  easely  explained  by  the  figures  and  that  summarize  the 
approach  we  are  taking.  George  is  in  Australia  and  the  secretaries  are  on 
strike.  I  had  to  cut  and  paste  because  1  cannot  find  the  originals.  I  hope  it 
will  serve  the  purpose. 


There  are  many  subtle  points  that  I  could  elaborate  on,  but  it  seemed  to 
me  that  you  would  prefer  more  pictures  than  text.  I  kept  the  text  to  a 
necessary  minimum  to  make  the  figures  understandable.  The  published  papers 
give  all  the  details,  so  it  would  be  unnecessary  duplication.  I  hope  the 
message  is  understandable. 


Thank  a  lot  for  your  help 


Manuel  Fiadeiro. 
Research  Scientist, 
Yale  University. 
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Ref .no.N00014-82-K-0371 

A)  THE  SEARCH  FOR  THE  LEVEL  OF  NO  MOTION 

The  Inverse  method  applied  to  the  mass  conservation  equations  Is  a 
powerful  method  to  obtain,  from  a  set  of  oceanographic  stations,  a  velocity 
field  in  geostrophic  equilibrium  obeying  the  conservation  constraints  imposed. 
The  solution  is  a  barotropic  field  that  represents  the  absolute  velocity  at  a 
given  reference  level.  If  there  are  less  constraints  than  station-pairs,  the 
problem  is  underdetermined  and  there  are  an  infinite  number  of  solutions.  In 
these  circunstances  one  has  to  choose  which  solution  to  adopt. 

It  has  been  the  practise,  for  the  purpose  of  calculating  the  absolute 
velocity  field,  to  assume  a  level  of  no  motion  within  the  water  column.  This 
heuristic  assumption  has  been  implemented  in  many  ways.  The  inverse  method 
provides  a  mathematical  formulation  and  a  quantative  way  to  assess  this 
assumption. 

If  a  level  of  no  motion  exists  and  the  mass  constraints  are  obeyed  it  is 
reasonable  to  assume  that  the  transport  imbalance  relative  to  that  level 
should  vanish.  Due  to  errors  in  the  data  we  can  expect  instead  a  small 
variance  from  zero.  We  proposed  a  procedure  in  which  the  reference  level  is 
changed  from  the  surface  to  the  bottom,  *o  search  for  the  minimum  magnitude  of 
the  transport  imbalance  vector.  This  procedure  gives  a  good  estimate  of  the 
mean  level  of  no  motion  and  in  most  cases,  acceptable  residues  for  the 
imbalances. 

We  can  also  use  the  projection  matrix  of  the  parameter  space  (the 
"resolution*'  matrix)  to  look  for  the  level  that  gives  the  smallest  of  the 
minimum  norm  solutions.  These  two  procedures  give  a  "best"  reference  level  no 
more  than  100m  apart.  In  most  cases  the  residual  error  and  the  variance  of  the 
solution  are  within  acceptable  limits  and  no  correction  is  necessary.  Examples 
are  shown  for  the  Coral  Sea  (Gascoyne  2/60,  1962)  and  Bermuda  triangle 
(Atlantis  215,  1955)datasets. 

The  Bermuda  triangle  needs  special  mention.  Here  we  divided  the  water 
mass  in  five  layers.  The  upper  two  layers  contained  shallow  stations  at  which 


Me  assumed  zero  velocity  at  the  bottom  (OV  stations).  Some  of  these  stations 
are  across  the  Gulf  Stream  were  this  assumption  Is  most  likely  to  be  false. 
Thus,  we  could  expect  large  Imbalances  for  those  two  layers.  We  then  calcu¬ 
lated  the  residual  transport  for  the  bottom  three  layers  Including  only  those 
stations  with  a  maximum  depth  near  or  below  the  reference  level.  The  search 
procedure  gave  us  a  best  reference  level  at  3500m.  The  imbalance  of  the  upper 
two  layers  was  then  attributed  to  the  0V  stations  and  an  inverse  solution  was 
found  for  them.  The  resulting  Gulf  Stream  transport  across  the  Ft.  Pierce 
section  became  comparable  with  the  measured  flow. 

As  a  final  note  we  restate  that  these  results  are  not  a  proof  that  a 
level  of  no  motion  exists.  They  only  show  that  there  are  solutions  compatible 
with  a  level  of  no  motion.  While  the  problem  is  underdetermined  there  will  be 
an  infinite  number  of  solutions.  The  inverse  solutions  are  projections  of  the 
true  solution  into  a  subspace  of  the  parameter  space.  We  have  access  only  to  a 
limited  representation  of  the  solution,  but  at  least,  this  representation 
should  vanish  if  the  assumption  is  true.  So  far,  we  have  always  found  a 
solution  compatible  with  a  level  of  no  motion  and  have  presented  the 
baroclinic  field  relative  to  that  level  as  the  most  probable  absolute  velocity 
field. 

Publications: 

On  the  determination  of  absolute  velocities  in  the  ocean  by  Manuel  E. 

Fiadeiro  and  George  Veronis.  J.  Mar., Res.  40(sup. ) ,159-182,1982. 

see  also:  Comments  on  "On  the  determination  of  absolute  velocities  in  the 

ocean"  by  L.L.  Fu  and  Reply  to  L.L.Fu  J.  Mar.  Res. ,42,259-262,1984. 

Circulation  and  heat  flux  in  the  Bermuda  Triangle  by  M.E.  Fiadeiro  and 

George  Veronis.  J.  Phys.  Ocean. ,13,1158-1169,1983. 
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B)  THE  INVERSE  PROBLEM  WITH  TRACER  MIXING 


Figure  1  shows  the($  and  PV  (f  |£>)  vertical  sections  along  the  meridian 
50°W,  Atlantis  cruise  ATL  229.  To  the  left,  St. 5416  Is  off  the  Grand  Banks  and 
to  the  right,  St. 5471  Is  off  the  French  Guiana.  The  section  closes  a  large 
section  of  the  Atlantic  Ocean,  the  Caribbean  and  the  Gulf  of  Mexico. 

The  assumption  that  properties  are  conserved  by  layers  cannot  hold  for 

such  large  bodies  of  water.  To  the  south  (right),  the  water  masses  present  a 

characteristic  combination  of  G\  and  PV  found  nowhere  else.  To  the  north,  the 

0  8 

deep  water  between  1000  and  3000m  has  an  homogeneous  PV*10  .  Between  Sts. 

5436  and  5449  a  minimum  PV  at  $3<40.7  marks  the  body  of  18°  water  that  has 
annual  contact  with  the  surface  further  east  and  to  the  north.  In  such 
circunstances  one  would  be  hard  pressed  to  choose  conservative  layers.  It  is 
obvious  the  different  water  masses  are  being  mixed  inside  the  region  closed  by 
the  section. 

This  example  points  out  the  need  to  include  mixing  effects  in  the  formu¬ 
lation  of  the  inverse  problem  if  we  hope  to  obtain  any  meaningful  results.  To 
that  purpose  we  undertook  a  theoretical  study  of  the  mixing  effects  on  the 
inverse  procedures. 

We  started  our  study  with  two-dimensional  tracer  distributions  generated 
by  a  model.  We  used  an  extremely  simple  flow  field  (a  constant  flow  from  left 
to  right)  and  a  constant  eddy  diffusion  coefficient.  The  direct  problem  was 
solved  in  a  grid  of  16x64  points.  For  the  inverse  problem  a  small  region  of 
6x6  points  is  sampled  and  treated  as  observed  data. 

In  figures  2  to  7  the  original  flow  field  would  be  represented  by 
equidistant  streamlines  with  a  value  of  0  at  the  bottom  and  4  at  the  top.  We 
show  the  result  of  different  procedures. 

Fig. 2  shows  solutions  of  underdetermined  systems,  using  the  equations  for 
only  one  tracer;  Cj  and  C£  are  different  tracers  and  (5,7)  and  (23,6)  refer  to 
the  sampling  region  of  the  direct  solution.  The  solutions  are  entirely 
different  from  the  input  flow  field.  The  patterns  change  with  the  tracer  used 
and  region  sampled.  These  are  minimum  norm  solutions;  we  can  retrieve  the 
original  flow  field  by  using  two  tracers  or  one  tracer  and  additional  correct 
constraints  to  make  the  problem  well  or  overdetermined. 

Fig. 3  represents  an  Interesting  case.  We  added  Incorrect  information 
(u«G)  to  the  system;  a)  with  a  weight  of  10”*,  b)  with  a  weight  of  10”*.  The 


streamlines  change  with  the  weight  of  the  velocity  constraints.  As  the  weight 
decreases  the  pattern  tends  to  the  minimum  norm  solution  shown  In  fig. 2c.  This 
is  in  contrast  with  the  cases  where  we  added  correct  Information;  the  correct 
streamfunctlon  was  retrieved  Independently  of  the  weight  of  the  velocity 
constraints.  This  can  prove  to  be  a  very  useful  technique  In  practise: 
“orthogonal11  Information  compatible  with  the  tracer  Is  accepted  by  the  system 
but  Incorrect  Information  that  conflicts  with  the  tracer  is  rejected  by  the 
system. 

In  practise,  data  have  errors  and  the  parameters  are  mean  properties  of 
the  field  for  the  scale  of  the  observations.  We  added  random  noise  and  sampled 
the  model  data  at  every  other  grid  point.  In  these  circunstances  the  model 
parameters  give  residuals  to  the  equations  and  the  inverse  procedure  finds  a 
solution  with  even  smaller  residues.  The  question  is  to  determine  how  far  is 
the  solution  from  the  expected  parameters  of  the  model.  The  underdetermined 
cases  were  hopeless,  the  residues  can  be  made  to  vanish  but  the  solution  has 
no  physical  significance.  Fig. 4  shows  solutions  of  strict  well -determined 
systems.  We  added  to  the  tracer  equations  the  condition  v=0  at  the  bottom  and 
top  of  the  sampling  region.  Fig. 5  shows  solutions  when  we  impose  v*0  all  over 
the  sampling  region  (the  system  is  overdetermined).  In  fig. 6  we  show  solutions 
with  a  no  vorticitty  condition  for  the  9  interior  points  (^£'0»  and 

fig. 7  solutions  with  two  tracers.  As  can  be  seen,  under  certain  conditions, 
overdetermined  system  (least-squares  solutions)  can  give  reasonable 
approximations  of  the  expected  flow  field.  As  a  rule  of  thumb,  we  found  out 
that  if  the  inverse  of  the  condition  number  of  the  coefficient  matrix  is 
smaller  than  the  mean  residue  square  of  the  equations,  the  solutions  are  very 
good  approximations  of  the  expected  velocity  field;  otherwise,  one  should  be 
extremely  careful  to  interpret  to  results. 

Regarding  the  application  to  oceanographic  data,  we  can  know  very 
precisely  the  condition  number  of  the  inversion  matrix,  once  the  problem  is 
set  up;  but  we  lack  a  good  estimation  of  the  errors  in  the  data.  We  are 
presently  working  on  this  problem. 


Unique  combinations  of  ^  and  PV  indicate  mixing  effects  in  the 
portion  of  the  ocean  closed  by  the  section. 
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Obtaining  velocities  from  tracer  distributions  by  M.E.  Fladelro  and 
George  Veronls.  submitted  to  0.  Phys.  Ocean. 

Inverse  methods  for  ocean  circulation  by  George  Veronls.  Advanced  Study 
Institute  Lectures,  Scripps  Institution  of  Oceanography.  In  press 
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